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SUMMARY AND INTRODUCTION 


This report contains results of a follow-on effort on the NASA Advanced 
Transport Operating Systems (ATOPS) Technology Studies contract (NAS1-1 6202). 
This previous task, sponsored jointly by NASA and the FAA, was a study of 
curved, descending approaches using the microwave landing system (Ref. 1). 

The major goal of the present study is to generalize the lateral and vertical 
guidance algorithms to be applicable to any desired approach path defined by 
waypoints specified in a general form, rather than being restricted to the 
particular approach paths previously studied. 

Waypoint data for the Washington National river approach was used with the 
generalized guidance laws in an MD-80 simulation program. Results from this 
program consist of performance data, which included lateral and vertical 
tracking errors, bank and pitch angle time histories, and ground track 
trajectories. The river approach was defined by seven waypoints and contained 
a pitchover to a 3 deg glideslope during the first turn. 

Ten cases were simulated for the river approach. Four cases started at entry 
to the approach and included two turns and the pitchover. These cases were 
run to study the effects of MLS noise, steady winds, and a speed change. The 
remaining six cases contained four turns at a constant glideslope down to 100 
feet altitude. These runs provided data on the effects of MLS noise, steady 
winds, turbulence, windshear, and a light-weight configuration. 
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CONCLUSIONS AND RECOMMENDATIONS 


The goals of the present study (to generalize previously developed guidance 
algorithms and to generate performance data for a particular simulated 
approach) have been achieved. A generalized waypoint database was established 
and a method developed that allowed almost any approach path to be specified. 
It became necessary to completely reformulate the lateral guidance algorithm 
to accommodate this new path definition concept. Only minor changes to the 
vertical algorithm were required. These revised guidance laws do work well 
with this new path definition. Tracking errors for ten different cases were 
generated and were within acceptable limits. 

There are two major areas that have been identified for potential follow-on 
tasks. They are missed approach/precision departures and fixed-base simulator 
evaluations. For the departure studies, the following specific tasks are 
recommended: 

Develop path definition for departures using only the back azimuth 
angle and range information. 

Establish lateral and vertical guidance algorithms that can be used 
with both missed approaches and precision departures. 

Simulate the missed approach/precision departures and evaluate system 
tracking performance. 

For the simulator evaluations, preliminary analytical efforts indicated the 
following tasks: 

Formulate studies to establish system architecture with regard to a 
fixed-base simulator. 

Coordinate MLS/guidance computer hardware/software interfaces. 

Develop a guidance computer in connection with real time simulator 
implementation. 
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MLS OPERATIONS 


Basically, the MLS consists of ground stations, an airborne receiver, and a 
display /control panel for the pilot. A type of operation, called an ILS 
"look-alike", is where conventional approaches would be flown using the MLS 
equipment. Initially, existing aircraft can be reconfigured to operate in 
this mode since it requires minimum airborne equipment changes. To realize 
the additional benefits of MLS over ILS, an on-board guidance computer is 
added to process steering signals for complex, curved paths. A database of 
waypoints (distances, altitudes, speeds) is also required as part of the 
system. The guidance computer uses the MLS signals and waypoint data to 
generate steering commands for landing and departures. These commands can be 
used to drive a flight director for pilot-in-the-loop operation. Alternately, 
the signals can drive the autopilot directly to provide autoflight operation. 
This latter mode is the one for which performance data were generated in the 
present study. 

COVERAGE VOLUME 


Figure 1 is a sketch of the existing ILS coverage characterized by its narrow 
beam. Signals for landing (provided by localizer and gli deslope transmitters) 
are available only within the indicated volume. In contrast. Figure 2 
represents the MLS coverage showing the possibility of curved paths, segmented 
glideslope, path lengthening, as well as ILS look-alike approaches. This 
coverage is provided by an azimuth transmitter (AZ), an elevation transmitter 
(EL), and precision distance measuring equipment (DME/P). Both the AZ and EL 
antennas radiate scanned beams to provide angular coverage in the indicated 
volume. The DME/P (colocated with the AZ equipment) is omnidirectional and 
provides precision range information to the aircraft. 

Also shown in Figure 2 is the coverage for takeoff and go around. This 
coverage is provided by a back AZ antenna and the same DME/P as used for 
landing. Elevation angle information is not available for takeoff or go 
around. 
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APPROACH GEOMETRY 


Angle and range variables for landing are defined in Figure 3. The AZ-DME/P 
equipment is located 1000 feet from the end of the runway (a 10,000 feet 
runway was assumed for this study). The EL transmitter is actually located 
400 feet from the runway centerline, but to simplify the simulation it was 
assumed to be at the coordinate origin. 


The position of the aircraft MLS antenna relative to the ground equipment can 
be calculated from > e EL » and R DME * For the assumed rectangular 
coordinate system, the following equations result: 


= - r dme sin 9 


AZ 


x ' = d cos e c . 
m EL 


cos 2 0 el (Rq ME 


d2 sin2 0 EL } ‘ y m 2 


z ; = - 


x' 2 + y' 2 
m 


tan 0 


EL 


( 1 ) 

( 2 ) 

(3) 


In Eqn(2), the term d is the distance from the ground coordinate center to 
the AZ/DME location. Altitude to the CG is positive for the geometry shown, 

so h'=-z' . In a later section, these position coordinates will be 

mm 

corrected for the offset of the aircraft MLS antenna from the CG. 
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APPROACH PATH DESCRIPTIONS 


WASHINGTON NATIONAL RIVER APPROACH 


Six waypoints, and the glide slope intercept point, were defined for the 
Washington National river approach. These waypoints are identified by name in 
Figure 4. One additional waypoint was added so the first part of the 
simulated flight would be at a constant altitude. Pitch over to a 3 deg 
glideslope is specified at ENTRY. Coordinates for the waypoints are given in 
Table 1. The radii shown were used in the simulation, although other values 
could be chosen for the same waypoint sequence. 

The numbered waypoints in Figure 4 are the points of turn entry and turn exit 
computed by the algorithms for circular transitions. Waypoint 1 was selected 
at an altitude of 100 feet, waypoints 2 and 3 define the end and beginning of 
the turn at FAF18, etc. Waypoint 14 was arbitrarily selected as the starting 
point for the simulated runs. These 14 waypoints are used by the lateral and 
vertical guidance laws to define the desired approach path. The definition of 
these waypoints and their numerical values will be discussed in a later 
section. 
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GUIDANCE CONCEPTS 


Although the present study uses the same basic MD-80 simulation program as 
reported in Reference 1, the lateral and vertical subroutines were completely 
rewritten. This revision was due to the generalization of the previous 
lateral and vertical guidance laws. A maximum of 20 waypoints can now be 
accommodated and the approach can be simulated. Flow charts and equations for 
these guidance concepts are discussed in detail in the following sections. 

WAYPOINT TRANSFORMATION 


Consider waypoints that are defined as the end points of straight-line 
segments or legs. If a radius is chosen for the circular part of the 
transition between waypoints, then the beginning and end of the turn (for a 
constant speed approach) are defined. These two points are uniquely 
determined by knowing three adjacent waypoints. Figure 5 depicts a general 
three-point geometry. The coordinates shown are in the ground plane, where an 
x-y coordinate pair is defined for each waypoint. The numbering of waypoints 
was chosen in decreasing order as the approach progresses. Waypoints 
WP ( i s ) , WP 1 ( i s -1 ) , and WP ' ( i $ -2 ) are the three input points. The 
integer 1 identifies the particular set of points and is used in the 
transformation calculations. The points WP ( i s ~l ) and WP( i s ~2) are the 
ones to be determined. Three angles Ap A 2 , and A^ of Figure 5 are 
easily determined from the input points. 


A-j = tan 


-i y'ty 1 ) - y(y 

x' (i $ -l ) - x(i s ) 


rad 


^2 ~ tan 


-i y'(y 2 > - y(yn 

x 1 ( i s ~2 ) - x ( i S ~1 ) 


rad 


A3 = it + I R (i s -l) (A-|-A 2 ) rad 


where: I^i^-l) = +1 for a right hand turn 

= -1 for a left hand turn 


(4) 


( 5 ) 
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The desired points on the circular arc that are tangent to the adjacent legs 
are calculated from: 


x(i $ -l) = x'(i $ -l) - R(i S ~1 ) cos A^tan (A 3 /2) (6) 

y(i $ -l) = y'(i s -l) - R(i $ -1) sin A^tan (A3/2) (7) 

x(i s -2) = x‘(i s -l) + R(1 S -1) cos A 2 /tan (Ag/2) (8) 

y(1 s -2) = y'(1 s -l) + R(i s -D sin A 2 /tan (A3/2) (9) 


The above process is repeated for each set of three adjacent waypoints until 
the transformation is complete. Results for the river approach are shown in 
Table 2. Once the coordinates in the x-y plane are known, the length of each 
leg is calculated. Altitudes at each waypoint are computed by plotting 
altitude vs the length of each leg. For the river approach, this plot is 
simply a continuous line at the 3 deg gli deslope. Altitudes shown in Table 2 
were computed in this manner. 

LATERAL GUIDANCE EQUATIONS 


The lateral guidance geometry was established by assuming the approach 
sequence contained a linear leg followed by a circular leg. A general two-leg 
segment is shown in Figure 6. Each leg is identified by the waypoint 
designation. Leg i lies between WP(i ) and WP(i -1). The integer i 

3 5 5 5 

starts at its maximum value (14 for the river approach) and is decreased by 2 
until i $ =2. Legs 2 and 1 are both straight-line segments (between WP(2) and 
GPI18 in Figure 4). The flow chart in Figure 7 shows the general program flow 
for the lateral guidance subroutine. Reference is made to the equations that 
are programmed in the subroutine as explained in the following sections. 

Along-Track Distance Computations. 

Each leg distance is computed by the lateral guidance subroutine, which also 
calculates the distance to go along the curved path (measured in the x-y 
plane). This distance is computed as a function of time. Centers of the 
circular arcs are also computed for use in the roll steering command. 
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( 10 ) 


First, the desired course angle, (i $ ) 

-i - y( 1 s ) 


i|> D (i s ) = tan' 


x ( i 1 ) - x ( i s ) 


, is calculated: 

rad 


where 0 < ^ D (i s ) < 2-rr 

Next, the linear leg distance is computed, followed by the total chord length 
for the entire circular leg. The arc length for the circular leg is computed 
from this chord length: 


L (y = [( x( v i} ■ x(i s } ) 2 + ( y(i s- 1) ' y{i S } )" 

T (i S ~1 ) =|^x(i s -2) - x(i $ -l)j 2 + ^y(i s -2) - y(i s -l)) 

-1 / * 5c t ( 1 s“ 1 ) \ 

D C ( V 1) - 2RC1 S .1) .in 


1% 


( 11 ) 

( 12 ) 

(13) 


Arc center for the curve starting at WP( i g -l ) is found from: 

X C ( i S~l ) = x(i S -1) " W 1 * R (V 1} sin W (14) 

y c (i s -D = yO s -U + I R( i s- 1 ) R(i s _1) cos W {15) 


Leg distances for legs i s and i $ -l (including i $ =2) are computed from: 


d lg^s^ MV 
d lg^s - 1 ^ = D c (l s _1) 


(16) 

(17) 


Total distance to go is found by summing each leg distance for 1 < 2 . 


VZ[ d l(V + W>] (i s <2 > 

n s 


The final two linear leg distances are added In separated. 

4 


% 


D c (l) = [ x(l ) 2 + y(l) 2 
°LG {2) = D L (2) 

d lg (D = d c (D 


r n =T n + 0,(2) + D c (i) 
°TG °TG ,L L 


(i s = 2) 


(18) 
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A quantity, D^, that will be used in computing the distance to go as a 
function of time is initialized at this point, and is initially set 


D 






(19) 

(20) 


Lateral Complementary Filter. 


Lateral position coordinates from the MLS, equations (1) and (2), are 
processed using complementary filters. These filters are the same as those 
used in Reference 1. For completeness, the equations for the filtered 

X A 

position coordinates ( ,y' ) and the filtered rate coordinates ( x y ) 

are repeated below: 


e 


x 


e 


y 


£ x 


A 



A 



•A. 



A 



X* " X m 

m m 


K 


y 1 

- y' 

1 

■ 

AT 

X 

m 

1 

y 1 

AT 


x' + AT I 0.654e 
m- I > 


y ;_ + AT 1 0.654e x 


x m + AT 

m- 


[ 0.129.x 


y m . + AT] 0.129e v 


< 500 ft 

< 500 ft 


♦M 

+ 0. 125e- + V 

X 

+ 0.1 25e^ + y 


(21) 

( 22 ) 

(23) 

(24) 
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In the above equations, the minus subscript denotes the first past value of 
the variable. 


The ground track angle, vp G , is computed from the rate estimates. The 
position estimates are corrected for the MLS antenna being displaced ^ feet 
from the aircraft CG longitudinally and feet vertically: 


i|> G * taif’(y m /x m ) 0 < i|) G < 360° (25) 

x m = - *, cos I’G (26) 

\ = h - \ s,n ♦g < 27 > 

h m = h ' - i sin 6 + £ cos 0 (28) 

m m i 2 


Ground speed and the ai stance to the next waypoint are calculated from: 



Linear Leg Steering Commands. 

Path deviation normal to the linear leg is calculated, and this error signal 
is used in forming the roll steering command. In Figure 8, the aircraft 

A /S. 

position { * m ,y m ) is shown with a tracking error Ay . The following 
equations define the roll steering command, 4>sTR * used to maintain track 
along the linear leg. 

° 2 = I x m ■ x( v ] 2 + l\i - *<y ] 2 

Ay = [ y m - cos W'> -[ x m ■ x(1 S‘' ) 1 sin W') < 31 ) 

Dj =[[> 2 - Ay 2 ]' 2 (32) 
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** = *G - W 


Ay = Vg sin A\p ‘ 

(33) 

^STR = " G Ll (Ay + G L2 A '^ 

(34) 

where = 0.0275 and G L2 = 18.18 


°TG = D 1 " d i 

(35) 


Turn Anticipation. 

Logic is used to anticipate the turn before the aircraft gets to the next 
waypoint. This action was found necessary to reduce the tracking error during 
transition to the circular leg. A distance, D^, is computed based on a four 
second lead time. As the aircraft nears the next waypoint, the logic switches 
the roll steering command to a ramp to start the turn. This event occurs at 
point 1 in Figure 9. When the aircraft passes the line at point 2 in the 
figure, the roll steering command switches to the circular tracking signal. 

As this switching occurs, the variable Dj is reset to the remaining distance 
to go. Also, the quantity Dj is set to zero at this time. These variables 
are now reinitialized properly for use in the next set of calculations. 

Circular Leg Steering Commands. 

Figure 10 shows the tracking error, AR , during the turn. The following 
equations define the roll steering command used for this leg. The reference 
course angle, , defines the desired path direction during the turn. It is 
used to calculate the deviation from the ground track angle, Eqn (35), which 
in turn defines the error rate term AR • 


C 


2 


AR 


( X <V 1 > - x m)) 2 + ( y( 1 S- 1 > - y m ) 2 


I R (i s -l) R (i s -1) 
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j, = tan' 1 — - 1s 11 , 0 < * R < 2i> 

' y c (i S- 1) " 

AR = V G sin (4» 6 - ^ R ) 


^STR tan 


-1 


32.2 R(i $ -1 ) 


- G^gAR - G^AR 


where G^g =0.01 and G^ =0.1 


(37) 

(38) 


(39) 


The distance along the circular leg is computed next. This distance, Dj, is 
measured from the last waypoint to the point projected along the desired 
circular arc (Figure 10). The total distance to go is then computed from this 
leg distance. 



I R <V 

-1) + AR 

2 + R 2 (i $ -1) - C 2 

2 

hv 1 ) 

+ AR 

R(i s -D 


D I ° ROs-UA 

°TG = D 1 ■ d i 


rad 


(40) 

(41) 


Turn Termination. 


As the aircraft approaches the next waypoint, logic is used to lead the 
transition to the linear leg similar to the turn anticipation. Point 1 in 
Figure 11 is calculated using a four second lead time. At this point, a 
signal is computed based on the next (linear) leg. This signal has the same 
form as the steering command, Eqn (33). The circular steering signal, 

Eqn (37), is still being used to command the bank angle during this time. At 
point 2, the linear signal and circular steering signals are equal in 
magnitude. At this time, the steering command is switched to the linear one, 
Eqn (33). Abrupt changes in the command signal are avoided by this logic. 
When point 3 is reached, the appropriate variables are reinitialized for the 
linear leg, and steering continues using Eqns (31) -(34). 
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VERTICAL GUIDANCE EQUATIONS 


Generalization of the vertical guidance algorithm discussed in Reference 1 was 
made during this study. The concept employed allows the vertical tracking 
error to be computed by knowing only the parameters that define the current 
leg in which the aircraft is flying. A constant gli deslope can be defined for 
each leg, or one gli deslope change can be executed per leg. 

One way to visualize this concept is to consider "straightening out" the 
approach path to lie in a single plane. The altitude is displayed vs the 
distance to go as depicted in Figure 12. A typical leg is shown as identified 
by waypoints WP(&) and WP(£-1). Altitudes at each waypoint are known, and the 
distance Dj is computed from Eqn (32) or Eqn (40) depending on whether a 
linear or circular leg is the current leg. Each waypoint has associated with 
it two glideslopes. (They may have equal or different values.) Consider the 
situation in Figure 12 where the aircraft is below the first glideslope for 
leg a . The desired altitude, h d , for this glideslope is computed from 

h D = h(A) - Dj tan G $ U,1) (42) 

After the aircraft has progressed to the region of the second glideslope as 
shown in Figure 13, the desired altitude is given by 

h D = h(£-l) + [ D l U) - Dj ] tan G $ U,2) 

Deviation from the desired altitude is defined by 

Ah ' h D - h m 

The error rate signal is derived from rate and acceleration signals. 

Figure 14 is the block diagram that contains a complementary filter that 
provides this error signal. Filtered altitude rate, hp , is a signal 
available in the MD-80 autopilot. It is used to develop filtered altitude 
rate error, ^ , as defined by the following difference equations that are 

equivalent to the conti nous system of Figure 14: 


(43) 

(44) 
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h = - [ h p - h F _ ] /AT 

/\ 

Az-| = C-| Az-|_ + Cgh 

Az£ = C-j AZ2_ + C2 1 Ah - Ah_ 

/N 

Ah = AZ-j + AZ2 


C 

C 

C 


1 

2 

3 


= e~ A "^ T (t = 4.0 sec) 

1 -AT/x 
= t 1 - e 


= 1 


g-AT/x 


/AT 


( 45 ) 


(46) 


As shown in Figure 14, the pitch command signal to the autopilot, 6^ , is 
composed of three components. The pitch steering signal is a combination of 

the altitude error and rate signals: 

/\ 

^STR = G V2 A ^ + G V4^ (47) 


where Gy 2 = 0.0588 and Gy^ = 0.294 


An integral term, 0 CI , is used on the Ah signal, and a predict term, 0 p , is 
used when switching from one glideslope to another. This predict term is 
computed from the vertical flight path angle, y , at the time of switching. 

0 CI = 9 CI- + G V3 AhAT (48) 

where Gy^ = 0.00294 



0 p = 6p. + T - G s (£,2) 
0 CE = 0 STR + 6 CI + 6 P 


(49) 

(50) 

(51) 


Switching to the second glideslope is determined by logic based on the second 
glideslope parameters. This signal is computed while the aircraft is still 
tracking the first glideslope. When this signal changes sign, switching 
occurs with a smooth transition to the second glideslope. The flow diagram 
for the vertical guidance subroutine is shown in Figure 15. 
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SPEED CONTROL 


Programmed speed and flap changes were simulated during leg 10 of the river 
approach. This choice was somewhat arbitrary, but it was made so that the 
speed change would be on the linear leg after the second turn. The result is 
the speed transient settles out before the next turn. 

Speed control is simply implemented by commanding a step flap change to the 
flap servo, and a step speed change to the autothrottle. All existing logic 
and limits in the MD-80 autothrottle remain unaltered. The result is a flap 
change from 15 deg to the desired setting of 40 deg at the max servo rate. A 
speed change from 165 knots to the target speed of 140 knots accompanies this 
flap change. 
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SIMULATION DESCRIPTION 


MD-80 SIMULATION PROGRAM 


A modified version of the MU-80 batch simulation program was used for the 
river approach studies, where the ILS subroutines were replaced with coding 
that simulated the MLS angle and range signals. Routines for lateral and 
vertical guidance and speed change were developed and added to the program to 
provide a closed loop simulation capability. Figure 16 shows the relationship 
of the basic elements of the simulation. The roll and pitch inner loops of 
the MD-80 autopilot accept commands from the MLS guidance computer. Surface 
commands generated within the autopilot provide the desired lateral and 
vertical path tracking control for the approach. 

Another mode of operation is possible where the autopilot is replaced with a 
flight director. A pilot would then fly the MLS approach and would make the 
necessary lateral, vertical, and speed maneuvers. A simulator is planned to 
be used in future pilot-in-the-loop evaluations of MLS approach procedures, 
pilot displays, and path tracking performance. The last sections of this 
report discuss some aspects of hardware interfaces and software requirements 
of the fixed-base simulator. 

WIND AND TURBULENCE MODELS 


Included in the simulation are options for winds, whose magnitudes are a 
function of altitude, and for turbulence as additive terms to the winds. The 
factor that is a function of altitude is denoted W and is defined by 

W = 0.43 log h + 0.35 (52) 


where h is the altitude in feet. Input values for S , S , and S ( the 

a y z 

wind magnitudes in knots) when the airplane is on the ground are used to 
compute the three wind components at altitude: 


W x - W S x 

W y = « s y 


(53) 
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Turbulence is a random variable that may be added to the above winds. White 

noise is passed through first order filters whose parameters are functions of 

the wind components as shown in Figure 17. Outputs of the filters are denoted 

U_, V. and W„ and are used to form the total simulated winds: 

9 9 9 

M Tx * W x + U g 

w Ty ‘ w , + v g 

Wj z W Z + Wg 

Aircraft performance characteristics limit wind components at touchdown to 
maximum values of 25 knot headwinds, 10 knot tailwinds, and 15 knot crosswinds. 


W1NDSHEAR 


One case was simulated in which a severe windshear was present during the 
final approach at Washington. This shear was based on the Logan/ Iberia 
incident of 1973. The data (published by Aviation Week in 1975 (Reference 2)) 
was modified for simulation purposes and appears in Table 3. Linear 
interpolation was used to define the shear between the altitudes listed. 

MLS NOISE MODELS 


Angle and range signals as output from the MLS receiver will contain noise. 
For simulation purposes, additive noise has been included in the azimuth, 
elevation, and DME range variablese AZ ,e EL ,&R DME respectively. The variables, 
as used in Equation (1) - (3), were modified to include this additive noise: 

0 AZ = e AZ(TRUE) + 0 AZ(NOISE) 

0 EL = e EL(TRUE) + 6 EL(N0ISE) 

R DME = R DME(TRUE) + R DME(N0ISE) 

The International Civil Aviation Organization (ICAO) has specified noise 
levels for which the MLS must meet. These levels, as defined by 2 a or 95% 
probability values, are given in Table 4. Values for the low frequency path 
following error (PFE) and the high frequency control motion noise (CMN) are 
specified for AZ, EL, and DME/P. 
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The ICAO- specified noise levels, if actually present in the MLS signals, would 
result in unacceptable aircraft activity when coupled to flight control 
systems capable of accurate path control. Flight measurements indicate that 
the noise is less. Consequently, assumed levels of reduced value were 
simulated. These reduced levels were based on actual flight test data 
performed in France in 1983 (Reference 3). Results of these tests were used 
to define the gains and frequencies shown in Figure 18. For simulation 
purposes, ia or 68 %, probability values are used. The reduced values shown in 
Table 4 were used to obtain the noise model of Figure 18. Azimuth and 
elevation levels are significantly lower than the ICAO values, whereas the DME 
level remains unchanged. 

During the course of the present study, two baseline cases were run with no 
noise for comparison purposes. All other cases contained the reduced noise. 
For certification purposes, the ICAO values must be used. However, it was 
felt that noise models, based on actual flight test data, would yield results 
representing levels encountered in actual practice. 
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SIMULATION RESULTS 


Two sets of runs were executed for the river approach. In the first set, four 
runs simulated entry to the approach at approximately 3000 feet altitude. In 
the second set, six runs simulated the final approach starting at 1700 feet 
and ending at 100 feet. Pitchover and a speed change were programmed at the 
higher altitude to study their effects on tracking. Tracking accuracy was of 
particular interest for the second set of runs in the presence of turbulence 
and shear. 

ENTRY APPROACH RESULTS 


Results for the initial path approach are presented for four cases as defined 
in Table 5. Specifically, the path starts at WP(14) (Figure 4) and ends on 
leg 10 after WP(10). Leg 14 is at a constant altitude. Pitchover to a 3 
degree glides! ope occurs at ENTRY. Thereafter, a constant 3 degree descent is 
maintained. 

Case 1. 

This is a baseline case with no noise or other disturbance. Bank angle time 
history as shown in Figure 19 has the waypoints identified along the time 
axis. Note that the bank angle starts to change before WP( 1 3) and WP (11). 

Turn anticipation accounts for this effect since the waypoints are defined at 
the time of switching from the linear steering command to the circular 
steering command. 

Pitchover occurs at 20 sec, with the pitch response shown in Figure 20. There 
are pitch corrections at each waypoint to account for the turning aircraft. 

As expected, a smooth transition between waypoints occurs, accompanied by a 
smooth pitchover to the desired gli deslope. This part of the approach poses 
no particular problems to the guidance algorithms. 

The lateral tracking error is shown in Figure 21. Discontinuities in this 
trace are due to switching of the steering signal at each waypoint. This plot 
is a composite of the Ay error, Eqn (31), and the AR error, Eqn (36). These 
errors, and their respective rates, are used to define the lateral steering 
command. Since this case has no disturbances, these errors are inherent in 
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the MU-80 and are due to the aircraft dynamics, autopilot inner loops, and 
limits in the autopilot. 

Figure 22 shows a time history of the vertical tracking error. At pitchover, 
the aircraft is about 65 feet below the desired altitude. This transient 
settles out to within a few feet for the remainder of the run, with some 
correction during the turns. Figure 23 is a plot of vertical descent rate. 
After transients, this rate is close to the -12.4 feet/sec nominal value for 
this descent. 

A computer-generated plot of the ground track for Case 1 is shown in Figure 
24. The dashed line connects the named waypoints and the solid line is the 
ground track. At the scale shown, the tracking errors are not obvious. This 
type of plot will be useful for cases with wind and turbulence when the track 
is plotted using an expanded scale. 

Case 2. 

Results for the MLS noise case differ from the baseline primarily in the 
vertical plane. Lateral tracking is minimally affected. (See Figures 25 and 
26.) Pitch response (Figure 27) shows the effect of the vertical tracking 
error (Figure 28) which contains the noise components. It should be noted 
that drawing conclusions using a single noise run is not valid. The trend, 
however, in pitch response after the turns does show approximately a +0.5 
degree dispersion about the nominal. Vertical tracking error shows a _+20 foot 
dispersion during the same period. Vertical descent rate. Figure 29, has a 
+2 ft/sec spread about nominal after the pitchover to 3 degrees. 

Case 3. 

A headwind and a crosswind were simulated as a function of altitude for this 
case. The magnitudes at threshold were 25 knots for the headwind and 15 knots 
for the crosswind. At the initial altitude of 3146 feet, these values (as 
modified by Eqn (52)) become 46 and 28 knots, respectively. 

Tne bank angle time history of Figure 30 shows a smaller angle achieved during 
the first turn (between WP(13) and WP(12)). This reduction is a result of the 
decreased ground speed due to the wind. About a 302i reduction in this speed 
had occurred at the end of the first turn as compared to the no-wind case. 
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Without wind, a 50 foot error to the left of the desired track developed at 
the end of the first turn (Figure 26). 

Since the bank angle was smaller with these particular winds, the tracking 
error was actually less (15 feet) because the aircraft did not turn as sharply 
(Figure 31). Had the wind direction been reversed, the opposite situation 
would have occurred. Tracking errors are very dependent on the relative wind 
direction, and in some cases the errors are reduced while in other cases they 
are increased. At the end of the second turn (WP(10)), the wind case shows a 
smaller error than the no-wind case. A 40 degree course change had occurred 
between the first and second turns. Consequently, the wind direction was 
favorable to cause the aircraft to track more closely. 

This particular wind condition was actually selected for the final approach 
cases. There is a restricted area near WP{5) and a crosswind from the right 
is a worst-case condition. 

Case 4. 

This last case for the initial approaches is a programmed speed change and a 
configuration change. Initially, the speed is 165 knots and the flap setting 
is 15 degrees. On leg 10, the speed is reduced to 140 knots and the flaps are 
extended to 40 degrees. The major effect on bank angle is to increase the 
commanded angle because of the higher speed (Figure 32). Both turns occur at 
165 knots. The lateral tracking error (Figure 33) is generally smaller as 
compared to the constant speed case, whereas the vertical tracking error is 
about the same. The speed and flap change commands occur slightly before 80 
seconds, and the speed change is smooth (Figure 34). 

TERMINAL AREA APPROACH RESULTS 


Results for the remaining six cases as defined in Table 5 are presented. For 
these cases, the initial position of the aircraft is on leg 10 and descending 
along a 3 degree glideslope. Four turns are executed during this final 
approach. Two additional effects studied for these cases are the addition of 
turbulence to the steady winds and a windshear profile. 
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These six cases were simulated down to a 100-foot decision height. At this 
point, the tracking errors were recorded and the speed error noted. Results 
are shown in Table 6. It must be pointed out that results for the cases 
containing noise or turbulence are random in nature. To obtain an accurate 
picture, one should process many runs to obtain a statistical dispersion. 

Time did not permit this processing during the present study. An alternate 
way of assessing performance of random disturbances is to compare time 
histories of the various cases. Relative dispersions over a reasonable time 
period can give some insight into the tracking performance. 

Case 5. 

This case is similar to Case 1 in that no noise or disturbances are present. 
Bank angle time history is shown in Figure 35, where the waypoints are 
identified on the time axis. Note that waypoints 5 and 4 are spaced only a 
few seconds apart. Consequently, the turns during this time are almost an 
S-turn type of maneuver. 

Pitch response. Figure 36, shows a +1 degree variation in angle due to the 
coupling from roll. For the initial approach cases, the variation was about 
+0.5 degree. More closely-spaced turns cause more pitch activity. 

Lateral tracking error, shown in Figure 37, is on the same order of magnitude 
as Case 1. Vertical errors of Figure 38 reflect the larger pitch variation 
and are about twice those in Case 1. Due to the expanded scale of the descent 
rate (Figure 39), the excursions appear to be large. Actually, they are on 
the same order of magnitude as the descent of Case 1 after the pitchover has 
occurred. The ground track for this case from the initial point to threshold 
is shown in Figure 40. 

Case 6. 

As in Case 1, the effect of MLS noise is seen in increased vertical errors 
with about the same lateral errors (Figures 41-45). Dispersion in pitch 
tracking error is about double that of the no-noise case (+15 ft vs +7 ft). 
This same level of MLS noise is used in the remaining cases. 
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Cases 7 and 8. 


The same headwind and crosswind (scaled by altitude) as used in Case 3 was 
used in Case 7. In addition, turbulence, as defined in Figure 17, was added 
to the wind for Case 8. 

Lateral tracking errors are compared for these two cases in Figure 46. 

Although the shapes of these two curves are different than the no-wind case 
(Case 5), the excursions are about the same. For the vertical plane. Cases 7 
and 8 are compared in Figure 47. Here, the comparison to the no-wind case 
shows an increase in excursion of a factor of 3 to 4. The steady wind values 
used in the simulation were somewhat severe. Since the RMS levels of the 
turbulence are functions of these steady wind levels, turbulence is also 
severe. Nevertheless, the expanded scale ground track shown in Figure 48 
shows a reasonable trace. Case 8 results produced no problems in flying near 
a restricted area, which is located about 1500 feet from leg 4 (the last leg 
before FAF18). 

Case 9. 

Windshear components used for this case are shown in Table 3. Positive values 
represent tailwinds and crosswinds from the right. The main effect of this 
windshear is to appreciably increase the bank angle required for the turns (up 
to 25 degrees as shown in Figure 49). The lateral tracking error, Figure 50, 
is about the same as Case 6. Vertical tracking error (Figure 51) increased to 
about 25 feet during the first part of the approach, with an increase up to 30 
feet near the change to a headwind (between 300 and 400 feet altitude). This 
increase represents a factor of 2 as compared to the error in Case 6. 

Case 10. 

This case represents a minimum-weight landing configuration. It is the same 
as Case 6 except for the weight and speed differences (See Table 5). As noted 
in Table 6, the error at the 100 foot altitude point are the same as the 
full -weight case. Peak error (Figure 52) was slightly higher (200 feet) at 80 
seconds as compared to a 50 feet peak for Case 6. 
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FIXED BASE SIMULATOR CONSIDERATIONS 


A logical extension to the current MLS guidance system simulation studies is a 
man-in-the-loop, real-time simulation evaluation. This task would involve 
taking the algorithms developed during the present study and revising them as 
necessary for use in a real-time environment. A fixed-base simulator that 
includes MD-80 autopilot and Performance Management System (PMS) hardware is 
available for this purpose. Interfaces between this hardware and the computer 
simulation program need to be established. 

The first step in identifying these interfaces is to note what inputs and 
outputs are associated with the algorithms. Table 7 is a signal computation 
diagram containing a review of the major input/output and computations 
performed by the lateral computer program. Table 8 is a similar diagram for 
the vertical guidance program. These tables indicate the basic signal 
generation and inputs necessary for their computation. 

Figure 53 is a simplified block diagram showing the interconnections between 
the MLS guidance computer and the other aircraft systems. The new MLS 
guidance functions will have to interface properly with the hardware that is 
included in the simulator and with the existing software. Tasks were 
identified in the Recommendations section that relate to this interfacing 
follow-on effort. 
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TABLE 1 - WASHINGTON NATIONAL RIVER WAYPOINTS 


x y Alt Turn Radius 


WAYPOINT NAME 

Ft. 

Ft. 

Ft. 

Ft. 

ENTRY 

-44603. 

36527. 

3146. 

9000 

CABIN 

-40951 . 

26491 . 

3146. 

8000 

WINDY 

-1 8426. 

15014. 

1353. 

7500 

KEYBR 

-16425. 

6990. 

948. 

6500 

MEMBR 

-10827. 

4014. 

628. 

7000, 

FAF18 

-5500. 

0 . 

285. 

7000. 

GPI1 8 

0 . 

0 . 

0 . 

_ _ _ 
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TABLE 2 - TRANSFORMED RIVER WAYPOINTS 


WAYPOINT NUMBER x(i s ) y(i g ) h{i s ) 

Ft. Ft. Ft. 


1 

-1 908. 1 

0 . 

100. 

2 

-3157.9 

0 . 

165.5 

3 

-7370.5 

1409.5 

402.39 

4 

-1 0387. 

3682.4 

600.33 

5 

-1 1 31 4. 

4272.7 

657.98 

6 

-1 3870. 

5631 . 5 

809.68 

7 

-17125. 

9798. 1 

1095.0 

8 

-17599. 

11698. 

1197.7 

9 

-21471. 

16566. 

1533.8 

10 

-38143. 

25060. 

2514.4 

11 

-42029. 

29453. 

2829.1 

12 

-43632. 

33860. 

3074.8 

13 

-46928. 

38155. 

3146.0 

14 

-48535. 

39280. 

31 46. 0 
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TABLE 3 - VALUES FOR SIMULATED WINDSHEAR 


ATTITUDE TAILWIND/HEADWIND CROSSWIND 


FT. 

KNOTS 

KNOTS 

5000. 00 

23.00 

-26.00 

1000.00 

23.00 

-26.00 

800. 00 

22.00 

-25.00 

600.00 

20.00 

-24.00 

500. 00 

18.00 

-23.00 

400.00 

12.00 

-17.00 

300. 00 

6.00 

-12.00 

200.00 

-3.00 

-4.00 

100.00 

-6.00 

-2.00 

0.00 

-4.00 

-2.00 
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TABLE 4 - MLS NOISE MODELS ( 2 a or 95% PROBABILITY) 


ICAO 

VALUES 


REDUCED 

VALUES 


AZIMUTH 


ELEVATION 


DME 


PFE = 0.104 DEG 
(+20 FT) 

CMN = 0.0548 DEG 
(+10.5 FT) 


PRE = 0.133 DEG PFE = 100 FT 

(+2 FT) 

CMN = 0.0663 DEG CMN = 60 FT 

(+1.0 FT) 


PFE = 0.0312 DEG 
(+6 FT) 

CMN = 0.0208 DEG 
(+4 FT) 


PFE = 0.018 DEG PFE = 100 FT 

(+0.30 FT) 

CMN = 0.012 DEG CMN = 60 FT 

(+0.20 FT) 


PFE = PATH FOLLOWING ERROR 
CMN = CONTROL MOTION NOISE 
RUNWAY LENGTH = 10,000 FT 
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TABLE 5 - SIMULATED CASE CONFIGURATIONS 



SPEED 

WEIGHT 

MLS 

STEADY 


WIND 

CASE 

(KNOTS) 

(LB) 

NOISE 

WIND 

TURB 

SHEAR 

1 

140 

132,000 

NO 

NO 

NO 

NO 

2 

II 

II 

YES 

II 

II 

II 

3 

II 

II 

II 

YES 

II 

II 

4 

165/140 

II 

II 

NO 

II 

II 

5 

140 

II 

NO 

11 

II 

II 

6 

140 

II 

YES 

II 

II 

II 

7 

If 

II 

II 

YES 

II 

II 

8 

II 

II 

II 

II 

YES 

II 

9 

II 

II 

II 

NO 

NO 

YES 

0 

116 

90,000 

II 

II 

II 

NO 


PATH 


INITIAL 


FINAL 
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TABLE 6 - TRACKING ERRORS AT 100 FT DECISION HEIGHT 



VERTICAL 

LATERAL 

SPEED 


ERROR 

ERROR 

ERROR 

CASE 

FT 

FT 

KT 

5 

2 

50 

< 0.5 

6 

2 

50 

<0.5 

7 

12 

30 

3 

8 

12 

10 

9 

9 

20 

30 

3 

10 

2 

50 

<0.5 
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TABLE 7 - LATERAL GUIDANCE SIGNAL COMPUTATION DIAGRAM 


INPUTS 


OUTPUTS 


MLS Angles, Range: ®AZ’^EL’^DME “ ^ 


MLS Position Coordinates: Xm’ y m ,h m 


Waypoint Data: x, y 

Turn Radius: R 


Desired Course Angle: 

Leg Distance: \g 
T otal Distance To Go: ^qlq 


Inertial Acceleration 

From Existing ^ 

Filters: x,y 


Filtered MLS Position and 

/\ /\ /\ 


Date Signals: (x^y^) ,(x m 



Ground Track Angle: ^G 

/\ A ^ 

Transformed MLS Position: x m ,y (tl ,h lll 
Ground Speed: Vq 
C ourse Deviation and Rate: Ay, Ay 
Roll Steering Command: <f> STR 
Distance to Go: d tg 
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TABLE 8 


VERTICAL GUIDANCE SIGNAL COMPUTATION DIAGRAM 


INPUTS 


Waypoint Data: h,G $ I 

Leg Distance: n ( 

LG 7 

Altitude Rate and 

Acceleration From 
Existing Filters: hp,h 


OUTPUTS 


Desired Altitude: h D 
Altitude Error: Ah 


A 

Altitude Error Rate: Ah 


♦ 

Pitch Steering Signal: e STR 
Vertical Flight Path 
Angle: y 

Pitch Steering Command: e CE 
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APPROACH 

DIRECTION 


XcOs-D.YcOs-l) 


Figure 5. Waypoint Transformation Geometry 



APPROACH 

DIRECTION 




INITIALIZE 

? 


DEFINE ALL CONSTANTS AND SET 
INITIAL CONDITIONS. COMPUTE: 
DESIRED COURSE ANGLE, EQUATION 
(10); LEG DISTANCES, EQUATION 
(16)-(17); AND THE TOTAL DISTANCE TO 
GO, EQUATION (20). 


RETURN 


COMPUTE: COMPLEMENTARY FILTER VARIABLES, 
EQUATION (21M24); GROUND TRACK ANGLE, EQUATION 
(25); CORRECTED AIRCRAFT POSITION FOR MLS 
ANTENNA LOCATION, EQUATION (26M28); GROUND 
SPEED, EQUATION (29); AND DISTANCE TO NEXT 
WAYPOINT, EQUATION (30). 



IS 

LEG <2 
? 


^ LINEAR n 
OR 

CIRCULAR 
v LEG ^ 
\ / 


CIRCULAR 


LINEAR 


COMPUTE: LINEAR TRACKING ERROR, EQUATION (31); 
ERROR RATE, EQUATION (33); LINEAR STEERING 
SIGNAL, EQUATION (34); AND THE DISTANCE TO GO, 
EQUATION (35): 




Figure 7. Lateral Guidance Flow Digram (Sheet 1) 
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ALTITUDE , ALTITUDE 





*P 



r 1=4.0 SEC 
G V2 = 0.0588 
G V3 - 0.00294 
G V4 = 0.294 


Figure 14. Block Diagram of the Vertical Steering Algorithm 
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ENTER 



COMPUTE: COMPLEMENTARY FILTER VARIABLE, EQUATION 
(45); VERTICAL FLIGHT PATH ANGLE, EQUATION (48) 



Figure 15. Vertical Guidance Flow Diagram (Sheet 1) 
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AZIMUTH AND PITCH AND ROLL 



Figure 16. Block Diagram of MLS Pitch and Roll Guidance Systems 


WHITE NOISE 

O— 


ff WN ~ 10 


a V 2r 
1 + rS 


TURBULENCE, KT 

► 



FILTER PARAMETERS 



LONGITUDINAL 

LATERAL 

VERTICAL 

7<SEC) 

600/V tas 

600/V tas 

30/v tas 

<*KT) 

0.15 IW X I 

0.15 IWyl 

1.5 


Figure 17. Turbulence Model Used in Simulation 



SHAPING FILTER RECEIVER 

FILTER 


AZIMUTH 

u, = 0.001 RAD/SEC 
u h = 0.16 RAD/SEC 
a WN = 0.064 
(Trec = 0.02 DEG 


ELEVATION 
= 0.001 RAD/SEC 
a> h - 0.34 RAD/SEC 
^WN "* 0.024 
<Treq * 0.0097 DEG 


DME 

= 0.001 RAD/SEC 
w h = 0.245 RAD/SEC 

"WN “ 1547 
a REC = 53 4 FT 


Figure 18. MLS Noise Model Used in Simulation ( la or 68-Percent Probability) 
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PITCH ANGLE (DEG) , , BANK ANGLE (DEG) 



0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0 110.0 

TIME (SECONDS) 

Figure 19. Bank Angle, Case 1 (Initial Approach, Baseline) 



TIME (SECONDS) 

Figure 20. Pitch Angle, Case 1 (Initial Approach, Baseline) 
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LATERAL TRACKING ERROR (FT) 




TIME (SECONDS) 

Figure 22. Vertical Tracking Error, Case 1 (Initial Approach Baseline) 
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(id 000' I) BON 





VERTICAL TRACKING ERROR (FT) 


6.5 


o 

LU 
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O 

Z 
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X 

o 



0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0 110.0 

TIME (SECONDS) 

Figure 27. Pitch Angle, Case 2 (Initial Approach, MLS Noise) 



TIME (SECONDS) 


Figure 28. Vertical Tracking Error, Case 2 (Initial Approach, MLS Noise) 


5] 


^NGLE (DEG) 




IGLE (DEG) 


TIME (SECONDS) 

Figure 31. Lateral Tracking Error, Case 3 (Initial Approach, Wind) 
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Figure 33. Lateral Tracking Error, Case 4 (Initial Approach, Speed Change) 






Figure 36. Pitch Angle, Case 5 (Final Approach, Baseline) 
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VERTICAL TRACKING ERROR (FT) 
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Figure 39. Vertical Speed, Case 5 (Final Approach, Baseline) 


LRTERRL MLS GUIDRNCE 



DOWNRANGE (1,000 FT) 



Figure 40. Ground Track, Case 5 (Final Approach, Baseline) 
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Figure 50. Lateral Tracking Error, Case 9 (Final Approach, Wind Shear) 
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Figure 51. Vertical Tracking Error, Case 9 (Final Approach, Wind Shear) 
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Figure 53. Interface Between MLS Guidance Computer and Existing Aircraft Systems 
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